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Abstract 

Non-perturbative renormalization group approach to the dynamical chiral symmetry breaking is 
(N- an effective method which can accommodate beyond the ladder (mean filed) approximation. The 

■ usual method relying on the field operator expansion suffers explosive behaviors of the 4-fermi 

O , coupling constant, which prevent us from evaluating the physical quantities in the broken phase. 

Q \ In order to overcome this difficulty, we solve the fiow equation directly as a partial differential 

equation and calculate the dynamical mass and the chiral condensates. Also we formulate a beyond 
the ladder equation and it gives almost gauge independent results for the chiral condensates. 
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I. INTRODUCTION 



The Wilsonian renormalization group approaches to the continuum quantum field theory 



have been formulated 



-3j and developed in the past decades (see reviews [J-l7|). In these 
approaches, the renormalization group (RG) flow equation is defined by a functional differ- 
ential equation, the solution of which gives the partition function defined by the functional 
integral. Using this framework, we are able to obtain new approximation methods extracting 
the non-perturbative information of the partition function. Therefore we call this type of 
the RG "non-perturbative renormalization group" (NPRG). 

In order to solve the fiow equation approximately, we usually expand the equation in terms 
of the field operators and their derivatives. The derivative expansion has been applied to 
evaluation of the universal quantities such as critical exponents and anomalous dimensions. 

le expansion with respect to 



Although it is difficult 



For example, in the three-dimensional scalar field theories, t 
the field operator without the derivatives converges very well [8 
to confirm the convergence with respect to the order of the derivatives, the result of the 
expansion up to the 4-th derivatives agrees well with the Monte Carlo simulations 12|. 

The NPRG has been also applied to the analysis of the dynamical chiral symmetry 
breaking (D^SB) in the strong coupling gauge theories and its effective theory. As first 
noticed by Nambu and Jona-Lasinio ISj, the scalar 4-fermi operators become the source 
of the D^SB. From the viewpoint of the NPRG, when we lower the renormalization scale 
of the effective action, the strong gauge interactions induce the effective 4-fermi operators, 
which bring about the D^SB at the low energy scale [l^. Unless there are explicit symmetry 
breaking terms such as mass terms, the running 4-fermi coupling constant diverges at a low 
energy scale. This explosive behavior is nothing but a signal of the D^SB [4, 15|. 

On the other hand, the (improved) ladder Schwinger- Dyson (SD) equation has been 
frequently used for the analysis of the D^SB in the strong coupling gauge theories. This 
ladder approximation has the strong dependence on the gauge-fixing parameter. Moreover 
it is difficult to improve the ladder approximation systematically. However the framework 
of NPRG allows us to take account of the effects of the non-ladder diagrams in a systematic 
fashion Q, l^- 

Since the running 4-fermi coupling constants diverge at a critical scale as mentioned 
above, the RG flow can not go beyond the critical scale towards the infrared limit. In the 
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simple framework of NPRG which maintains the chirally symmetric structure of the effective 
action, we can not evaluate the physical quantities in the broken phase such as the dynamical 
mass of the fermion and the condensates of the fermion bilinear composite operator. 

To overcome this problem, we introduce the bare mass term, which also works as a 



source term of the chiral condensates 



18j. We may expect that the bare mass prevents 



the divergent behavior of the 4-fermi coupling constants and might allow us to effectively 
evaluate the order parameter of the D^SB at the infrared limit scale. The field operator 
expansion of the NPRG equation, however, does not converge well at least in the region 
of the bare mass as small as the current masses of up and down quarks. Consequently we 
will take another way of directly solving the NPRG flow equation as a partial differential 
equation (PDE) without relying on any field operator expansion. 

By the way, another method like the Hubbard-Stratonovich transformation has been 
used in many works to avoid the divergent behavior 19|-[2l|. In this method, the composite 
operators of fermions are partially transformed to scalar fields, one of which obtains the 
nonzero expectation value as a chiral order parameter. Introduction of those scalar fields has 
a merit that the meson physics can be argued simultaneously. However it is more complicated 
to evaluate the convergence of the physical quantities with respect to the operator expansion 
including the scalar fields. 

This article is organized as follows. In Sec. 2 we briefly review the flow equation for the 
effective average action. We introduce the basic truncation which projects the complete 
operator space onto the subspace relevant to the D^SB so as to solve the flow equation 
approximately. In Sec. 3 we examine the truncation method in detail and obtain two types 
of truncations: one corresponds to the ladder approximation, and the other contains the main 
parts of the non-ladder corrections. In Sec. 4 we explain how to evaluate the chiral order 
parameters. In Sec. 5, we solve the flow equation by using the field operator expansion, and 
examine the convergence of chiral order parameters. In Sec. 6, we directly solve the partial 
differential equation and show the results of physical quantities. In Sec. 7, we summarize 
our methods and results, and discuss further issues along the line of thought of this article. 
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II. NPRG FLOW EQUATION AND ITS APPLICATION 



A. Formulation 

In order to evaluate non-perturbative effects of the quantum field theory, we introduce 
the so-called "effective average action" jsl that interpolates between the bare action and the 
full quantum effective action. For this purpose, we define the generating function with the 
infrared cutoff supplied by the cutoff term A^Affi] as follows: 



Z^[J] := j p(]e-^'^--[^l-^^^[^]+^-^-^, (1) 



where VL represents various fields generically and the functional integral is regularized by 
the ultraviolet (UV) momentum cutoff Aq. The cutoff term AS'a is defined as the following 
mass term depending on the momentum: 

A5A[n] = hf{-v) ■ Ra{p) ■ m- (2) 

The regulator function Ra{p) suppresses the quantum fluctuations with the momentum 
lower than the infrared cutoff A. Therefore the regulator function satisfies 

hm i?A(p)>0, (3) 

and it implements the "coarse graining" in the Wilsonian method. 

The effective average action is defined by a slightly modified Legendre transformation: 

rA[$] = sup (^j ^-J- Wa[J]^ - A5a[$]. (4) 



It satisfies the following boundary conditions, 

Ta -> 5'bare, A Aq ^ OO 

Ta r, A -> 0, 

provided that the regulator function Ra{p) has the following properties, 

Ra{p) oo, a Aq ^ oo 
Ra{p) = 0, AVg2^0 



(5) 



(6) 
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The cutoff dependence of the effective average action is exactly given by the NPRG flow 
equation, 



dtT^m = -STr 



■dtRA 



(7) 



where we define the dimensionless scale parameter t = logAo/A. The right hand side should 
be called /3 functional, and consists of the second order functional derivative of the effective 
average action. 



^(2) 



{P,<l) 



(8) 



5$i(-p)(5$j(g)' 

It is considered as a matrix with respect to the momentum and the species of fields. This 
flow equation is firstly derived by Wetterich [Sj]. Because of the boundary condition ([S]), 
the flow equation interpolates between the bare action Sbare and the full quantum effective 
action F. 



B. Application to QCD 

Hereafter we consider the A^f-fiavor massless QCD with A^^-color. The bare action in 
Euclidean space is 

^bare = /{iC^r + ^ + ^9sm (9) 

where T" is the generator of the fundamental representation of SU(A^c)- This action has the 
chiral symmetry SU{N{)i^ x SU{Nf)ii, which is to be broken down to SU{N{)y dynamically 
by the strong gauge interactions. 

To evaluate the dynamical chiral symmetry breaking (D^SB) by using the non- 
perturbative renormalization group (NPRG) flow equation, we truncate field operators which 
are not essential to the D^SB from the complete operator space of the effective average ac- 
tion. Here we define the following truncated effective average action, 

^^[•^l = [{^^^'^^''^ + ^ + ^ ^^^^ + '3sA) ^ - V{^, V^; A) I , (10) 

where we use the co variant gauge with the gauge- fixing parameter ^, and do not represent 
the ghost sector for simplicity. This truncated subspace of the complete effective action is 
spanned by the operators of the bare QCD action and the fermion self-interaction operator 
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V{il!,'ijj; A), which we call the fermion potential. The discarded operators including higher 
derivatives or gluon fields may affect quantitative evaluation, but are not the essential sector 

le bare QCD and the 



19|. 



to drive the D^SB. Actually we will confirm that the operators of t^ 
fermion potential bring about D^SB using the NPRG fiow equation 

Here it should be noted that the fiow equation induces operators breaking the gauge sym- 
metry, such as the mass operator of the gauge boson, because the cutoff function exphcitly 
breaks the gauge symmetry. However we do not discuss this issue because we truncate out 
such operators. Eventually the flow equation of the gauge coupling constant agrees with that 
of the one-loop perturbation theory. We will concentrate on evaluating the flow equation 
for the fermion potential. 



C. Techniques to derive the NPRG flow equation 



In the rest of this section, we will explain some techniques to explicitly write down the 
NPRG flow equation ([7]) for the fermion potential. We denote the degrees of freedom of the 
fields as a vector $* = (AJ^, ■?/'*, ip). Then the dressed inverse propagator including the cutoff 
function is given by 

rf [<i>](p,g) + i?A(p,g) 
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rA|*] + ASA|*l 



\ Si>'(-p) I 

The regulator function is defined by 



Ra{p, q) 



\ 





















X 5ip-q), 



I 



(11) 



(12) 



where {D^^)f^{p) = " - |)) and5(p-g) = (27r)^5(^)(p - g). Functions, r{p) 

and r^(p) are defined to satisfy the properties ([3]) and ([6]). 

Next we explain how to calculate the "super trace" in the NPRG fiow equation ([7j). We 
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transform the flow equation as follows 



fl, 



dtT^m = dt isTr log 



ri'^ + Ra 



Here the symbol dt is defined by 

dt{ZAr{p)) 



+ 



dt{Z^r^{p)) 6 



(13) 



(14) 



(15) 



Za Sr{p) Z^ Sr^ip) 

where = J (f;;^- Then we split the inverse propagator matrix (ITT]) into submatrices as 
follows, 

'Mbb Mbf 
Mfb Mff 

where Mbb (Mff) corresponds to the second derivative with respect to bosonic (fermionic) 
fields, while Mbf and Mfb correspond to bosonic and fermionic field derivatives. Physi- 
cally the submatirx Mbb corresponds to the inverse free propagator of the gauge field, the 
submatrices Mbf and Mfb are the gauge interactions, and the submatrix Mff contains the 
inverse propagator of the fermion and the fermion self-interactions. Using this notation, we 
can rewrite the "super-trace log" in the flow equation (fT3|) into the following formula 22| : 

STr log M = -Tr log Mff + Tr log [Mbb - MbfMfp^Mfb] • (16) 

Note that another expression, 

STr log M = Tr log Mbb + Tr log [Mfb - MfbM^i^Mbf] , (17) 

has been often used in the NPRG analyses. Here Eq. (fT6|) is more appropriate for our 
purpose of improving the gauge parameter dependence. 

To derive the flow equation for the fermion potential, the matrix M in Eq. f ll5p should 
be evaluated by replacing the fields with their zero-momentum components. Replacing 
$(p) — > ^6{p), we have 

M(p,g) = rf [$](p,g) + i?A(p,g) 



[zFil + r)iD,') 



$(p)=$5{p) 

iZ^{l + r^)^ - 



dipdip 



iZ,^{l + r^)p - 



\ 



d^jjdtjj 



dip dip 



(18) 



X 6{p-q). 
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III. FLOW EQUATION FOR THE FERMION POTENTIAL 



A. Scalar 4-fermi operator 

The gauge interactions induce all possible fermion operators respecting the symmetry of 
QCD, which are enhanced by themselves when we lower the cutoff scale. Even in the trun- 
cated subspace of the effective action (ITU]) , we can not treat exactly all possible operators 
and interactions. Therefore, as an approximation of the flow equation for the fermion poten- 
tial, we project its full operator space onto a specific subspace and restrict the interactions 
so that we evaluate the D^SB most effectively. 



As in the QED with one flavor [16], the central operator for D^SB is undoubtedly the 
following scalar 4-fermi operator. 



where (/ = 1, ■ ■ ■ , A^f^ — 1) are the generators of the fundamental representation of 



model with A^f = 3, is invariant under the chiral transformation, SU{N{)-[^ x SU{N{)^. It is 
the only chiral invariant 4-fermi operator which gives corrections to the mass operator, and it 
becomes the relevant operator in the region of the strong gauge coupling constant. Therefore, 
for a flrst-step approximation, we project the operator space of the fermion potential onto 
the subspace spanned by polynomials in the scalar operator p. 

To project the flow equation onto the subspace deflned above, we determine the coefficient 
from all possible operators included in the full fermion potential. This is equivalent to count 
the coefficients of powers of {ipipY, even though this operator itself is not chirally invariant. 
It is due to the fact that {^ipY operator does not appear in chiral invariant operators other 
than powers of p. Eventually, we may work with a potential function in the simplest scalar 
operator a = ipip- 





■2 



(19) 



7=0 




V{ij,^P)->V{a). 



(20) 



Here we note that the original chiral symmetry is not maintained in this subspace, but the 
discrete chiral symmetry still remains: the Lagrangian (the fermion potential) is invariant 
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under the following discrete transformation, 

ip -)> 75V', -^75, o- -a. (21) 

The discrete chiral symmetry forbids the operators of the odd powers of a, such as a mass 
term. 

Next we pick up the interactions that are expected to be most important for the D^SB 
and for improvement of the gauge-fixing parameter (^) dependence. In Eq. (1T8|) . we further 
select the large- A^^c leading interactions in the fermion self-interactions, which leads to the 
following simplification: 

1 /2 

Applying the above approximation and the usual field renormalization, ip i^/^^ ^^^^ 

1/2 

— )■ A^j^/Z^ , to Eq. ( ITS]) and (ITB]) . we obtain the flow equation for the fermion potential, 
dtV{a; t) = -v^ad^V + jti 4 log S^\p) - ^ jti dt log [1 + A{p) + B{p)] , (23) 

where rj^ is the anomalous dimension of the fermion field. Functions A{p) and B{p) are 
matrices in the space of the color adjoint representation and the Euclidean Lorentz vector 
space as follows: 

[Ar'ip) = g'jT^^,S4p)^,D,,{p)T'ij, (24) 



{B)%{p) = gl,l;T'D,,{p)^,S^{-p)^,Ty. (25) 



Here we should pay attention to the difference between the two types of traces: tr acts the 
space of the fermion's representation, the Dirac spinor and the color fundamental represen- 
tation and the flavor fundamental representation, and tr acts the space of the matrix A, B. 
The propagators of the fermions and the gauge bosons including the cutoff function in Eq. 
flMj) and ([25]) are defined by 

^^^^^ - PM + d^v- ' ^ ^ 

^''"^^ p{p) y ^" ^ ^' p2 



(DTM = ^ V - (1 - 0^ , (27) 
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where P^pip) = (l + r^)^j9^, and P{p) = (l+r)p^. According to the field renormahzation, the 
anomalous dimension ?7^(= dtlogZ^) appears in Eq. fl25|) . and the gauge coupling constant 
is renormalized as Qs = ZiCjs/ Z^zlC^ . 

In the next two subsections, we will explain how to extract the scalar operators from the 
right hand side of Eq. fl23|l . 

B. Ladder approximation 

In Fig. [H the fiow equation ( 123|1 is diagrammatically expressed by using the "corrected 
vertex" defined in Fig. HJ We regard A and B as the ladder element and the crossed element 
of the diagrams, respectively. 




FIG. 1: Infinite series of diagrams expressing the flow equation (j23|) using the corrected vertex. 

^ = W ^ "XT 

\l V 

FIG. 2: The corrected vertex. The two diagrams in the right hand side correspond to A and B in 
Eq. (|23|) respectively. The curved arrows denote the direction of the momentum flow. 

The diagrams consisting only of the ladder element A are the ladder diagrams. Moreover 
we find that the diagrams consisting only of the crossed elements B are also the ladder 
diagrams if we rotate all the fermion lines of the diagrams to untangle the crossed gluon 
lines. Therefore the ladder approximation with only the ladder diagrams should be defined 
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by 



-trlog[l + A + B] = -J2 + BT 

n=l ^ 
ladder ly. Mr , 
n=l 



5^ ^^trA^ (28) 



n 

n=l 



To project Eq. (128|) onto the subspace of the scalar operators cr", we adopt the following 
approximation rule of picking up a, 

n 

i=l 

where we use T° = -^j^ Icoior- The right hand side of the above formula is nothing but 
the scalar part of the general Fierz transformation obtained by using the completeness of 
the space of the spinor and the color and the flavor. Thus the coefficient F„ in the above 
formula is given by 

Fn = tr [lspinorA°T°Oi ■ Ispinor A°r°02 " " " 1 spinor A °r°0„] . (30) 

According to this rule, the summation in Eq. fl25]) is calculated as follows: 
-1)" ~ 



CO 

^ n ^ 

n=l n=l 



n-i 



(31) 



where C2 is the second Casimir invariant of the SU (Nc) representation, C2 = Yla=i ^ T^-T"". 
Finally we obtain the ladder flow equation. 



dtV{a- 1) = -v^ad^V + jti dt log (^S;^\p) - C2 g, 



'Tmn • '''' 

where we rescaled V and a by common factor NcN{, V — A^^cA^f V and a — )■ AcA^f a. 
As for the regulator function, we adopt the following sharp regulator function, 

rsharp(p) = ZA-r{p) = Z^^- r^{p) = - 1- (33) 
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Performing the momentum integration in Eq. 
partial differential equation (PDE), 



, we obtain the ladder flow equation 



dtV{a; t) = -Tj^ad^V + — - In 



1 + A-M d^V + (3 + 



C2 tj 
4A2 



(34) 



Apart from the anomalous dimension term, this flow equation agrees with the local potential 
approximated Wegner-HougMon equatio,, and .t was proven that the flo^ equation gives 
the results equivalent to the improved ladder Schwinger- Dyson equation pj;]. Actually this 
is the reason why we call this approximated flow equation "the ladder" . 

Using the momentum scale expansion 23|, |2J], the anomalous dimension of the fermion 
field 1]^ = dt log is given by 



C29I 



A' 



3-e AW(t) 



7r2 V A2 + m2(t) 4 {A^ + 711^)^ 



where we define the running mass, 



m{t) = daV{a;t)\^=o- 



(35) 



(36) 



As will be seen in the numerical results in Sec. IVH the chiral order parameters given by 
the ladder flow equation (134|) strongly depend on the gauge-fixing parameter ^. 



C. Beyond the ladder approximation 

In order to improve the gauge dependence, we have to add the crossed element as well. We 
evaluate the flow equation fl23l) using the full corrected vertex in Fig. |2l which is calculated 
as follows: 

(A + B);l{p) = - p^l^ly^ ^T^T' (^ ^^^p^^^" (W757/.) + MD,^) ^ ^^^^ 

+ (term including [T^T^] ). 

Here we will ignore the commutator term for simplicity. Note that this ignorance is consistent 
with our approximation that we do not include diagrams with gluon self couplings. From 
the following discussion we expect that this ignorance does not induce the strong gauge 
dependence in the truncated subspace. 

We associate the corrected vertex with the gauge independent set of diagrams for the 
S-matrix in case of the Abelian gauge theory. In the non-Abelian gauge theory, the diagram 
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exchanging one gluon is also needed for the gauge independence of the S-matrix. On the 
other hand, such corrections can not be added directly due to the one-loop nature of the 
NPRG P function, and therefore the gauge dependence which would be canceled in the S- 
matrix appears in the commutator term of Eq. (137|) . In the NPRG, the correction by the 
diagrams exchanging one gluon is treated through the effective operators such as d^F^^'ip'-^yip. 
However such operators are not included in the truncated subspace. Thus we may omit the 
commutator term without suffering strong gauge dependence. 

Then, according to the general Fierz transformation fl291) . we can pick up the scalar 
operators as follows: 



^ ' -t"r(A + BY 2N^N{^ ' 



2n 



n 



a 



2P{P^ + M^)N,N,_ 



k=0 



n 
2k 



(2 + 4'=)M"-2'=P^ 



(38) 



Finally, adopting the sharp regulator function, we obtain the flow equation beyond the 
ladder approximation as the following partial differential equation: 



dtV{cr; t) = -rj^ad^V + — log 



4vr2 



+ ^log 



A2 + 52 



+ 



3A2G2 



A4 , 



1 + e 



d^VG 
A2 + d^V^ 



(39) 



where B = O^V + Csf^ and G = Caft^ 



IV. CHIRAL ORDER PARAMETERS 



Now we explain how to evaluate the two chiral order parameters, the dynamical mass of 
quarks and the chiral condensates (ipip), which are generated by the D^SB. 

In the framework of the NPRG, calculating the non-zero chiral order parameter is nontriv- 
ial because the NPRG flow equation maintains the chiral invariant structure of the effective 
action, which forbids appearance of the dynamical mass operator. ^ 

^ In the theories including the scalar fields whose symmetries spontaneously break, we can evaluate the 
nonzero expectation values of the scalar fields as the order parameters by searching the minimum point of 
its effective potential. In the theory we consider, this method cannot be used directly because the chiral 
order parameters are not such expectation values of the fields. Therefore, the scalar fields corresponding 
to the ^p'^p are introduced by methods like the Hubbard-Stratonovich transformation in many works. In 
this article, however, we adopt another method as will be seen in Sec. IVII 
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On the other hand the D^SB shows itself as a divergent behavior of the 4-fermi couphng 
constant, which is the source of the D^SB in the NJL model. We can define the /3 function for 
each operator by expanding the flow equation in powers of a. The /3 function for the 4-fermi 
coupling constant consists of itself and the gauge coupling constant, but does not include 
the higher dimensional operators due to the chiral invariance. Solving the RG equation, we 
obtain the flow of the 4-fermi coupling constant as follows: Lowering the cutoff scale A(t), 
the gauge interactions generate the 4-fermi operator, which enhances itself, and consequently 
the 4-fermi coupling constant diverges at a finite infrared scale A^. 

Because of this divergence, the RG flow can not go beyond the critical scale Ac toward 
the infrared limit in the chiral invariant operator space. Therefore the divergence seems to 
imply that the chiral invariant RG flow cannot exist at the cutoff scale lower than Ac, where 
the true RG flow might be in the chiral variant operator space including the mass operator. 
The relation between the divergence and the D^SB has been discussed in Refs. {4, 16, 25 1. 

Here, in order to go beyond the critical scale Ac and to effectively evaluate the chiral order 
parameters, we introduce the bare mass term, which explicitly breaks the chiral symmetry, 



in addition to the chiral invariant bare action: 



18| 



5'bare = 5'bare (invariant) - / rriQiljil). (40) 

J X 

When the cutoff scale A(t) lowers, the running mass m(mo; t), being mo at the initial scale 
t = 0, is rapidly enhanced around Ac by the 4-fermi interaction, and consequently the 
enhanced mass suppresses the (3 functions due to the decoupling effect. Therefore the 4- 
fermi coupling constant is expected to stay finite at the infrared scale. Taking the zero bare 
mass limit after solving the flow equation, the infrared limit mass becomes a chiral order 
parameter, the so-called dynamical mass of quarks: 

mdvn = lim limm(mo;t). (41) 

mo-5>+0 t-5>oo 

The bare mass mo also works as an external source for the composite operator ipijj, and 
its expectation value is evaluated as follows: 

(^^)^ lim lim ^^°i^°'^\ (42) 

mo-!>+oi-s>oo amo 

where Go(mo;t) denotes V{4',4'',t)\ii,=-4,=o- Here the infrared limit value (5o(mo; 00) corre- 
sponds to the Helmholtz free energy. 
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Prom the view point of the Helmholtz free energy, we may describe the relation between 
the DxSB and the divergent behavior of the 4-fermi couphng constant. Due to the chiral 
invariance of the theory, (5o(mo;t) is an even function of mo- Hence the derivative of 
Goimo'jt) at mo = 0, namely the chiral condensates, vanishes if it is an analytic function. 
On the other hand, the derivative of Go{fno; t) cannot be continuous at mo = 0, if the chiral 
condensates has the non-vanishing values. When the cutoff scale A(t) decreases, the initially 
analytic function Go(mo; t) changes into a non-analytic one at the critical scale Ap. Therefore 
its second derivative diverges at the scale Ap. Actually, in the large-N approximated NJL 
model, the 4-fermi coupling constant is almost equal to the second derivative, and therefore 
its divergent behavior shows the emergence of the D^SB at that scale Ac. 

V. FIELD OPERATOR EXPANSION AND ITS CONVERGENCE 

Here we attempt to solve the flow equation using the field operator expansion so that we 
evaluate the chiral order parameters introduced in the previous section. As seen in the end 
of this section, however, the field operator expansion does not work well in this model. 

In the subspace spanned by polynomials in cr, the (3 function for the 4-fermi coupling 
constant includes the 6-fermi coupling constant owing to the chiral symmetry breaking 
effect of the bare mass itlq. In general, the /? function for the 2n-fermi coupling constant 
includes the 2(n-|- l)-fermi coupling constant. Therefore we have to take into account of the 
infinite number of the 2n-fermi operators, and encounter the infinite tower of the coupled 
RG equations. 

Since such infinitely coupled equations can not be evaluated numerically, we stop the 
expansion at some maximum power N: 



where we call N the order of "truncation" . Then, expanding the flow equation in powers of 
0", we set up (A^+ l)-coupled RG equations where the (5 function for the 2A^-fermi operator 
does not include the 2 (A'"-!- l)-fermi coupling constant. The lager the truncation order N is, 
the better the solution of the truncated coupled RG equations approximates the solution of 
the original flow equation. Actually, this type of truncation approximation has worked well 
in many theories. 




(43) 



n=0 
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Let us see the truncation dependence of the dynamical mass mdyn calculated for the 
ladder flow equation f l34l) . In App. |X] we explain the input parameters and the running 
gauge coupling constant which obeys the one-loop perturbative /3 function. In Fig. |3], we 
plot the running mass m(mo; t) at the infrared limit ( t — > oo) for each bare mass, and show 
its dependence on the truncation order from = 2 to = 20. In the bare mass larger 
than about 0.02 GeV, the truncated solutions converge well with respect to truncation order 
A^. It also shows that the mass is dynamically generated by the D^SB since the infrared 
running mass is much larger than the corresponding bare mass. In the small bare mass 
region, however, the truncated solutions do not converge but diverge more badly for larger 
truncation order. Therefore we can not take the zero bare mass limit straightforwardly 
in this method with the field operator expansion. It is also difficult to make a reliable 
extrapolation. The non-ladder flow equation fl5^ shows the similar behaviors as the ladder 
flow equation. 
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FIG. 3: Dependence of the infrared limit of the running mass m(mo; t — > oo) on the bare mass mg 
and truncation order N . 

Here we show the results using another method of avoiding the explosive behavior of 4- 
fermi coupling constant. In Eq. f H3|l . the fermion potential is expanded around the vanishing 
value ((7 = 0). We expand it around the nonvanishing value do (7^ 0), 

^ 1 

y{^\t) = Y,-E^{k){a-a,r- (44) 

n=0 
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In this expansion, we can go beyond the critical scale without introducing the bare mass 
because the nonvanishing value ctq plays a similar role as the bare mass. Note that it does not 
work as an external source for the chiral condensates. The dynamical mass corresponding 
to Eq. (14T|) is given by the following limit, 



mdyn = lim lim Hi{ao;t). 

(70— >+0 t— ^OO 



(45) 



As seen in Fig. IU however, the expansion around the nonvanishing value does not converge 
in the small value of the expansion point (Tq. 




Iog(oo/On 



FIG. 4: Dependence of Hi{ao;t oo) on the expansion point do and truncation order A^, where 
the normalization scale of ctq is Cnorm. = 1 (GeV)'^. 



VI. SOLVING THE FLOW EQUATION AS PDE 

In Sec. IIVI and |Vl so as to evaluate the chiral order parameters, we have introduced 
two methods, introduction of the bare mass and expansion around the non-zero point, both 
of which are expected to avoid the divergent behavior of the 4-fermi coupling constant. 
Note that the expansion around non-zero point means that the fermion potential V{a] t) 
is evaluated at the non-zero point, a ^ 0. In these methods, however, the field operator 
expansion does not converge in the small value region of the bare mass or the expansion 
point, and we cannot obtain reliable chiral limit of the dynamical mass. 
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Obviously the poor convergence of the field operator expansion originates from the diver- 
gent behavior of the 4-fermi coupling constant, which corresponds to the second derivative 
of the fermion potential V{a;t). On the other hand, if the fiow equation is solved as a 
partial differential equation (PDE), the fermion potential is expected to behave analytically 
at least except for the origin. ^ Therefore, we will directly solve the flow equation as a PDE 
without the fleld operator expansion. 

In the practical calculation we solve the flow equation in terms of the mass function, 
M{a;t) = d(jV{a;t), because the numerical solution of the PDE for the mass function 
is more stable than that for the fermion potential. Moreover, for the numerical stability 
around o" = 0, a is transformed into the logarithmic variable, x = logcr/cTnorm where the 
normalization scale (Tnorm is set to be 1 (GeV)^. We do not directly treat the origin of a and 
consider the limit x — — oo of the mass function M{x; t) to calculate the dynamical mass 
as seen in Eq. 

Now we adopt the simple formulation of the grid method where the derivatives with 
respect to x are replaced with the flnite differences. The flnite differences are deflned by the 
7-point formula where for example the flnite difference corresponding to the flrst derivative 
dxM consists of values on 6 points. These points are nearest neighbors to the point where 
the derivative is evaluated. 

For the numerical calculation, we will set a flnite region for x, xi^ < x < x-^. In order to 
approximate the mass function at the end point xi^ to be the dynamical mass, we need to 
choose a small enough xl << —1. Near the boundaries, where we can not take full 6 points, 
we use the 5-point formula, the 3-point formula, and eventually at the very end we take only 
the next point. Usually these constraints enhance the numerical error near boundaries. It 
will be found that actually these low order approximated deflnition of the flnite difference 
near the boundaries does not induce the instability. In the practical calculation we choose 
the boundaries such that xl = —19 and xr = 0. 

Through the procedure explained above, the coupled ordinary differential equations for 
the discretized mass function on the grid are obtained, and we numerically solve the equa- 

^ Strictly speaking, such a function without the total analyticity cannot be a global solution of the PDE. 
Adopting the "weak" solution of the PDE, we can exactly solve the PDE and obtain physical quantities 
straightforwardlywithout any extrapolation. This method of the weak solution will be reported in a 
separate article |26| . 
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tions with respect to the dimensionless scale t using the 4-th order Runge-Kutta method. In 
Fig. Owe present the RG evolution of the mass function given by the ladder flow equation 
f l34|) . Here the Landau gauge ^ = is adopted, and the anomalous dimension is ignored, 
which is consistent with the local potential approximation (LPA). We can see that the mass 
function is dramatically increased from the vanishing value when lowering the cutoff scale 
A(t). Particularly the dynamical mass generation is observed below a critical scale, rather 
rapidly in a short range of scale t. The infrared limit is reliably evaluated, whose size reaches 
the Aqcd scale, and thus the D^SB occurs. If the chiral symmetry is not spontaneously 
broken, the mass function goes to be zero in the limit x — ?■ — oo. 




-16 -12 -8 -4 

log(f^/f^norm.) 



FIG. 5: RG evolution of the mass function M((t; t). The green dash lines denote the mass function 
at dimensionless scales t increasing by At = 0.3 step, and the red line denotes the mass function in 
the infrared limit, t oo. The mass function is dramatically increased from the vanishing value 
below a critical scale tc — 5.3. 

The mass function well converges to a certain value when x goes towards the end point 
xl • Therefore we can conclude that the approximated definition of the finite differences near 
the boundaries does not induce the instability of the mass function at xl, and the value can 
be identified with the dynamical mass mdyn- To evaluate the chiral condensates (ipip) by 
using Eq. fH2l) . we introduce the bare mass and calculate the free energy Gq through the 
evolution of V at the end point xl- As mentioned in Sec. IIIIBl the LPA ladder flow equation 



19 



with the Landau gauge has been proved to give the result equivalent to the improved ladder 
SD equation. Actually the two chiral order parameter obtained now agree well with the 



29| . which assures the total consistency of our 



ones obtained by the SD approach in Ref. 
method. 

We present the numerical results of the two chiral order parameters obtained from the 
two approximated flow equations, the ladder one fl34p and the non-ladder one (I39p . with the 
various values of gauge fixing parameter ^. Table [T] (Table [TTl) shows the numerical values 
of the dynamical mass (the chiral condensates) with or without the anomalous dimension. 



Here the chiral condensates are the renormalized ones at 1 GeV, {ipip)iGeV 29 1. 



TABLE L Results of the dynamical mass m-dyn. [GeV] . (L) ladder approximation; (NL) non- 
ladder approximation; (LA) ladder approximation with the anomalous dimension; (NLA) non- 
ladder approximation with the anomalous dimension. 



e (L) (NL) (LA) (NLA) 



0.888 


0.958 0.792 


0.860 


1 1.12 


1.14 


0.827 


0.860 


2 1.33 


1.30 


0.861 


0.844 


3 1.54 


1.44 


0.896 


0.821 


4 - 




0.931 


0.792 


5 - 




0.967 


0.757 



We show these results graphically in Fig. [6] and [71 Now we find that the gauge dependence 
of the results with the anomalous dimension (A. D.) are suppressed much better than ones 
without the A. D.. Besides, as for the results with the A. D., the gauge dependence of 
the chiral condensates obtained from the non-ladder flow equation is almost vanishing. It 
looks perfect and this is nothing but what we expected to get by adding the crossed gluon 
diagrams to go beyond the ladder. On the other hand, such great improvement by the non- 
ladder calculation is not seen for the dynamical mass. This is due to the fact that the chiral 
condensates are the on-shell quantities while the dynamical mass is not. Consequently there 
is no reason that the dynamical mass does not depend on the gauge fixing parameter. We 
confirmed that our non-ladder extended calculation respects the gauge independence almost 
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TABLE II: Results of the chiral condensates {{ipip)iGc'V /^f)^^^ [GeV]. The notation is the same as 
in Table HI 



C (L) (NL) (LA) (NLA) 

0.215 0.219 0.209 0.213 

1 0.267 0.264 0.221 0.219 

2 0.330 0.312 0.236 0.222 

3 0.408 0.364 0.255 0.224 

4 - - 0.278 0.225 

5 - - 0.304 0.225 



perfectly. 

Finally it should be noted that in the Landau gauge (^ = 0) the gauge dependent ladder 
result of the chiral condensates coincides with the almost gauge independent non-ladder 
extended one. This feature of the Landau gauge proves a folklore that the ladder approxi- 
mation is good particularly in the Landau gauge. 

~i 1 1 1 

a 

« ■ ■ ! 

ladder o 
ladder with A. D. 

non-ladder □ 

non-ladder with A. D. ■ 

J I I I 

1 2 3 4 5 

FIG. 6: Gauge dependence of the dynamical mass mdyn- 
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FIG. 7: Gauge dependence of the chiral condensates (C. C.) (V''i/')iGev/-^f • 
VII. SUMMARY AND DISCUSSION 

In this paper we have derived a new approximated flow equation beyond the ladder 
approximation so as to improve the dependence on the gauge fixing parameter. We have 
developed various methods to evaluate the two chiral order parameters, the dynamical mass 
of quarks and the chiral condensates. These methods are expected to avoid the explo- 
sive behavior of the 4-fermi coupling constant in the course of solving the flow equation. 
Within these methods, however, the field operator expansion, which is usually applied to 
solve the flow equations, shows poor convergence induced by the infrared singularity. Then, 
we stopped the field operator expansion, and solved the flow equation as a partial differ- 
ential equation by using the grid method. We have obtained the chiral order parameters 
successfully without any instability or extrapolation. 

As for the chiral condensates, we have seen that the gauge dependence of the non-ladder 
extended result almost disappeared. Therefore our non-ladder extended approximation al- 
most respect the gauge invariance. 

A next step for further improvement of approximation would be to take into account of 
higher order operators including the first derivatives. This improvement is implemented by 
replacing the coefficient of the kinetic term with the function of a, Z^{(j\t). Con- 

sequently we need to solve coupled partial differential equations in terms of the fermion 
potential V{a]t) and the kinetic function Z^{a]t). 
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Also we have seen that particularly in the Landau gauge the gauge dependent ladder 
result of the chiral condensates coincides with the almost gauge independent non-ladder 
extended one. This agreement might be related to a statement that only in the Landau 
gauge the ladder approximation satisfies the Ward-Takahashi identity. However, at finite 
temperature and chemical potential, this relation is broken. Therefore we would encounter 
the new phase structures by applying this non-ladder extended approximation to the hot 
and dense QCD in the framework of the non-perturbative renormalization group. 
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Appendix A: Running gauge coupling constant and input parameter 



In our truncated subspace, the /3 function for the gauge coupling constant, ctg = 91 /'^t^, 
agrees with the result of the 1-loop perturbation: 



/3o 2 

/vr 



(Al) 



where = 



|A^f. The solution of the RG equation is 



27r 



(A2) 



/3o i^qcd — t 

where the scale Aqcd = Aoe~*i'^'^ gives the scale of the infrared quantities such as the chiral 
order parameters. Here, in order to go beyond the scale lower than Aqcd, we introduce 



the infrared cutoff effect 



that the gauge coupling constant stop increasing at a proper 



infrared scale, which is naturally expected by the confinement. Adopting the cutoff scheme 



in Ref. 



291], the running gauge coupling constant is given as follows: 
^2ti 1 



/3o ^qcd — t ' 

27r 1 



TT 



2n 



qcd 



/3o (^ir — ti)(tqcd 



TT 



. tqcd — tir (tqcd — ^ir)^ 



t < tir 
ti < t 



(A3) 
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where we set a fixed dimensionless scale for ti to be tqcd + 1, and tir is left as an infrared 



cutoff scale parameter. Obeying Ref. 



291 ]. should be parametrized by Ajr as = tqcd — 



0.5 ■ (Air + 1), and we take the following parameter: 

Aqcd = 484 MeV, Ai^ = -0.5. (A4) 

By the analysis using the ladder Schwinger-Dyson equation, it is confirmed that the physical 
quantities such as the chiral condensates are not sensitive to the choice of the infrared cutoff 
scale parameter tij-. 

Finally, we discuss the initial condition of the fermion potential V{a;t). The initial 
cutoff scale Aq has to be large enough so that the obtained RG flow well approximates the 
renormalized trajectory starting from the ultraviolet limit, Aq — )■ oo. At the ultraviolet limit, 
the fermion potential vanishes because the effective average action agrees with the bare QCD 
action, or if it exists there it would be strongly suppressed soon by its higher dimensionality. 
In practical calculation, we take the vanishing fermion potential as the initial condition, and 
set the initial cutoff scale Aq to be large enough so that the infrared quantities does not 
depend on the Aq within a given numerical precision. Actually we set Aq to be the Z boson 
mass scale, Mz = 91.2 GeV. 
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